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Abstract 

Although rice (Oryza sativa) is the most Al-tolerant species among small-grain cereal crops, there is wide geno- 
typic variation in its tolerance to Al toxicity. A number of quantitative trait loci (QTLs) for Al tolerance have been 
detected, but the responsible genes have not been identified. By using chromosome segment substitution lines, this 
work found that Nratl, a gene encoding an Al transporter, is responsible for a QTL previously detected on chromo- 
some 2. Substitution of the chromosome segment containing Nratl from Koshihikari (Al-tolerant variety) by that from 
Kasalath (Al-sensitive variety) decreased Nratl expression and Al uptake and tolerance, but increased binding of Al 
to the cell wall. Nratl in Kasalath showed tissue localization similar to Koshihikari in the roots. Although Koshihikari 
and Kasalath differed in four amino acids in Nratl protein, Nratl from Kasalath also showed transport activity for Al. 
Analysis with site-directed mutagenesis revealed that these differences did not affect the Al-transport activity much. 
Furthermore, there was no correlation between Al tolerance and the open-reading-frame sequence differences in 
other rice varieties. On the other hand, there was good correlation between Nratl expression and Al tolerance; how- 
ever, sequence comparison of the promoter region up to 2.1 kb did not give a clear difference between the Al-tolerant 
and -sensitive varieties. Taken together, these results indicate that differential expression of Nratl is responsible 
for the QTL for Al tolerance on chromosome 2, although the mechanism controlling Nratl expression remains to be 
examined. 
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Introduction 

Aluminium (Al) toxicity is a major limiting factor of crop pro- 
duction on acid soils, but there is a wide difference in Al toler- 
ance between species and cultivars within a species (Ma et al, 
2002; Famoso et al, 2011). Among small-grain cereal crops, 
rice {Oryza sativa) shows the highest tolerance to Al toxicity, 
but genotypic differences also exist; usually Japonica varieties 
show higher Al tolerance than Indica varieties (Ma et al. , 2002). 
Furthermore, the relative degree of Al tolerance in the five sub- 
populations follows the order temperate japonica > tropical 
japonica > aromatic > indica = aus (Famoso et al, 2011). 

A number of quantitative trait loci (QTLs) for Al tol- 
erance have been identified by using different mapping 



Journal of 

Experimental 
Botany 

www.jxb.oxfordiournals.org 



populations and phenotyping methods (Ma and Furukawa, 
2002). Wu et al. (2000) identified four QTLs for Al tolerance 
on chromosomes 1, 3, 9, and 12, using relative root length 
as a physiological parameter in a random inbred mapping 
population derived from Azucena (Al-tolerant) and IR1552 
(Al-sensitive). Five and 10 QTLs, respectively, for Al toler- 
ance scattering on different chromosomes were detected in 
populations derived from Chiembau x Omon 269-65 and 
from CT9993 x IR62266 (Nguyen et al, 2001, 2003). On 
the other hand, Ma et al (2002) used relative root elonga- 
tion as a parameter and identified three QTLs for Al toler- 
ance on chromosomes 1, 2, and 6 in a population of 183 
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backcross inbred lines derived from a cross of Koshihikari 
(Al-tolerant) and Kasalath (Al-sensitive). Xue et al. (2007) 
mapped three QTLs for Al tolerance on chromosomes 1, 9, 
and 1 1 in a recombinant inbred line population derived from 
a cross between the tolerant japonica cultivar Asominori and 
the sensitive indica cultivar IR24 based on relative root elon- 
gation. Recently, a total of 48 distinct Al-tolerance genomic 
regions were detected by genome-wide association mapping 
based on relative root growth (Famoso et al. , 201 1). However, 
the genes responsible for these QTLs have not been identified. 

On the other hand, by using mutant approaches, a num- 
ber of Al-tolerance genes have been identified in Japonica 
varieties (Ma et al, 2014). ART1 (Al-tolerance transcrip- 
tion factor 1), a C2H2 zinc-finger type transcription factor, 
was reported to be involved in Al tolerance (Yamaji et al, 
2009). ART1 regulates at least 31 genes by binding to the 
core ra-acting element [GGN(T/g/a/C)V(C/A/g)S(C/G)] 
in the promoter of these genes (Tsutsui et al, 2011). The 
expression and localization of ART1 is not induced by Al, 
but the expression of downstream genes is upregulated by 
Al within hours. Seven ART 1 -regulated genes (STAR1, 
STAR2, Nratl, OsALSl, OsFRDL4, OsMGTl, and 
OsCDT3) have been functionally characterized. STAR1 and 
STAR2 encode a ATP-binding domain and a membrane- 
binding domain, respectively, of a bacterial type ABC trans- 
porter (Huang et al, 2009). The STAR1 STAR2 complex 
localized to vesicles transports UDP-glucose, which may 
be involved in cell-wall modification, resulting in decreased 
Al accumulation in the cell wall. OsFRDL4 is responsi- 
ble for the secretion of citrate in response to Al (Yokosho 
et al., 2011), while OsCDT3, encoding a small cysteine-rich 
peptide, shows binding activity with Al, thereby prevent- 
ing Al entering into the root cells (Xia et al, 2013). On the 
other hand, OsMGTl functions as a Mg transporter (Chen 
et al, 2012), and upregulation of OsMGTl could allevi- 
ate internal Al toxicity by enhancing Mg uptake. Nratl, 
a member of Nramp family, takes up trivalent Al at the 
plasma membrane (Xia et al, 2010), which is required for 
subsequent sequestration of Al into the vacuoles for final 
detoxification. Vacuolar sequestration of Al is mediated by 
OsALSl, a half-size ABC transporter localized to the tono- 
plast (Huang et al, 2012). Among these genes examined, 
OsFRDL4 showed a good correlation between its expres- 
sion level and Al tolerance (Yokosho et al., 2011), indi- 
cating that this gene may be responsible for the genotypic 
difference in Al tolerance. By contrast, there is no correla- 
tion between the genotypic variation of expression levels 
of STAR1, OsMGTl, OsCDT3, and OsALSl and Al toler- 
ance (Huang et al, 2009, 2012; Chen et al, 2012; Xia et al, 
2013), suggesting that they are involved in fundamental Al 
detoxification processes common in most rice varieties. 

Recently, genome-wide association detected a sin- 
gle, highly significant region on chromosome 2 that was 
unique to the aus subpopulation (Famoso et al, 2011). 
This region, which contained the Nratl candidate gene 
(Os02g0131800), is the same as the location of a QTL for Al 
tolerance on chromosome 2 detected previously in a popu- 
lation from Kasalath and Koshihikari (Ma et al., 2002). 



This QTL explained 7.3% of the variation of Al toler- 
ance. Physiological characterization showed that although 
Al-induced secretion of citrate from the roots was higher in 
Koshihikari than in Kasalath at higher Al concentrations, 
the difference between the two varieties was not significant 
(Mae? al, 2002). 

The present study investigated whether Nratl is responsible 
for the QTL detected by using chromosomal segment sub- 
stitution lines. Furthermore, expression level, tissue locali- 
zation, and Al transport activity were compared between 
Al-tolerant and -sensitive varieties. This work found that dif- 
ferential Nratl expression level is partially responsible for the 
genotypic difference in Al tolerance in rice. 

Materials and methods 

Plant materials and growth conditions 

Two chromosome segment substitution lines (SL204 and SL205) 
were provided by the Rice Genome Resource Center (http://www. 
rgrc.dna.affrc.go.jp/). In SL204, the segment from marker C1357 to 
G132 (0-60.3 cM of chromosome 2) containing Nratl was substi- 
tuted by the Kasalath segment in Koshihikari background, while in 
SL205, the segment from marker G132 to C747 (60.3-107.7 cM of 
chromosome 2) was substituted (Supplementary Fig. SI available 
at JXB online), which was used as a negative control. Rice seeds 
(Kasalath, Koshihikari, SL204, SL205) were soaked in tap water 
overnight at 30 °C in the dark and then transferred to a net floating 
on 0.5 mM CaCl 2 in a 1.5-1 plastic container. Seedlings were grown 
for 4-7 d at 25 °C. Similar size seedlings were selected and used for 
the following experiments. 



Evaluation of Al tolerance 

Six rice seedlings (5-d-old) per each genotype were exposed to 
0.5 mM CaCl 2 containing 0, 30, or 50 uM Al (pH 4.5) for 24 h. Root 
lengths were measured with a ruler before and after treatments. 
Relative root elongation was calculated as follows: (root elongation 
with Al) / (root elongation without Al) x 100. Six roots were meas- 
ured for each treatment. 



Al determination in cell sap and cell wall 

For determining Al accumulation in the root tips, 5-d-old seedlings 
(Kasalath, Koshihikari, SL204, SL205) were exposed to 50 uM Al 
(pH 4.5) for 8h and then root segments (0-1 cm, 20 roots each) were 
excised after washing three times with 0.5mM CaCl 2 . To obtain 
root cell sap, the root segments were put in ultra-free-MC centrifu- 
gal filter units (Millipore) and centrifuged at 3000 g for lOmin at 
4 °C to remove apoplastic solution. The roots were then frozen at 
-80 °C overnight. The root cell sap solution was obtained by thaw- 
ing the samples at room temperature and centrifuging at 20 400 g 
for lOmin. The residual cell wall were washed three times with 70% 
ethanol and then immersed in a 0.5ml of 2M HC1 for at least 24h 
with occasional vortex. The Al concentration in solution was deter- 
mined by inductively coupled plasma mass spectrometry (ICP-MS, 
7700X, Agilent Technologies). 



RNA isolation and gene expression analysis 

To examine the expression pattern of Nratl, rice seedlings (7-d- 
old) were exposed to different Al concentrations (0-50 uM) for 
8h. Root tips (0-1 cm) were excised. All samples with three repli- 
cates were subjected to RNA extraction. Total RNA was extracted 
using an RNeasy Mini Kit (Qiagen). Total RNA (1 |xg) was used for 
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first-strand cDNA synthesis using a Superscript II kit (Invitrogen), 
following the manufacturer's instructions, with an oligo(dT) 12 _i 8 
primer. Expression was determined with SYBR Premix Ex Taq 
(Takara) by Mastercycler ep realplex (Eppendorf). The primer 
sequences used for reverse-transcription PCR were as follows: 
Nratl, forward 5'-GAG GCCGTC TGCAGGAGAGG-3' and 
reverse 5'-GGAAGT ATCTGC A AGC AGCTCTG ATGC- 3 ' ; and 
ART1, forward 5'-CAGTGCTTCT CGTGGGTCTT-3' and reverse 
5'-CCTG TGCGTGAA GAACCACT-3'. HistoneHS (forward 
5'-AGTTT GGTCGCTCTCGATTTCG-3'; reverse 5'-TCAAC 
AAGTTGACCACGTCACG-3') was used as an internal control. 



Immunohistological staining 

Immunostaining with an antibody specific to Nratl was per- 
formed, obtained by immunizing rabbits with the synthetic peptide 
N-MEGTGEMREVGRETLHGG-C (positions 1-18 of Nratl) 
(Xia et al, 2010). After the roots of Koshihikari and Kasalath 
were exposed or not to 30 iiM Al for 12h, the root tips (2 mm) 
were excised. The procedures for immunostaining were the same as 
described previously (Yamaji and Ma, 2007). 

Yeast experiments 

The cDNA fragment containing the entire open reading frame 
(ORF) for Nratl derived from Koshihikari or Kasalath was ampli- 
fied by PCR using the primers 5'-GGTAC CAAAAT GGAAGGG 
ACTGGTGAGATGA-3' and 5'-CTACATGGAAGCATCGGCA 
A-3' (forward and reverse, respectively). The fragment was first 
cloned into the pGEM-T vector (Promega). After sequence confir- 
mation, Nratl cDNA was excised with Kpnl and Not! for cloning into 
pYES2 (Invitrogen). Site-directed mutants of Nratl (Koshihikari) 
were generated on the Nratl-pYES2 construct by PCR using the 
following synthetic oligonucleotide primers: common primers as for 
amplification; for mutation of Nratl (Koshihikari): E120K, forward 
5--GAGTACCCGCACTATGTCAACA-3' and reverse 5'-TGTTG 
ACATAGTGCGGGTACTC-3';V326I,forward5'-GATCTCCCCA 
AAATATTCCTGAGTAAAAGAGGTGTGG-3 and reverse 5'-AC 
TCAGGAATATTTTGGGGAGATCAAGCTCGGTG-3';T500M, 
forward 5'-CAGGAA GGACATGGTGGCGACGTACGTGC 
CG-3' and reverse 5'-ACGTCG CCAC CATG TCCTTCCTGA 
ACGTCAGGT-3'; and V5 1 5A, forward 5'-GTGG CGACG TACGT 
GCCGGTGC-3' and reverse 5'-GCACCGGCACGTACGTCGC 
CAC-3'. Using combinations of these primers for PCR, four mutant 
types of Nratl (Koshihikari) were generated. Mutation in each 
amplified fragment was confirmed by DNA sequencing. These frag- 
ments were then cloned into pYES2 (Invitrogen). The resulting plas- 
mids were introduced into Saccharomyces cerevisiae BY4741. 

Al sensitivity evaluation on agar and Al uptake in liquid culture 
were performed as described by Xia et al. (2010). For Al sensitivity 
evaluation, Nratl (Koshihikari), Nratl (Kasalath), and vector control 
pYES2 were transformed into S. cerevisiae BY4741 and then grown 
on solid media (LPM without uracil; MacDiarmid and Gardner, 
1996) containing 0, 200, or 300 uM Al buffered with 5mM succinic 
acid. For Al uptake in liquid culture, transformants were selected 
on LPM without uracil and grown in synthetic complete medium 
without uracil. Cells at midexponential phase were harvested and 
transferred to LPM without uracil. Cells were cultured for 2h, fol- 
lowed by addition of A1C1 3 to the medium to a final concentration of 
50 uM Al. After 4h of incubation with shaking, cells were harvested 
by centrifugation. The pellets were washed three times with milli-Q 
ultrapure water (Millipore) and digested with 2 M HC1. The concen- 
tration of Al in the digest solution was determined by ICP-MS. 

Western blot analysis of Nratl 

Total protein (5 iig) extracted from the yeast expressing different 
alleles of Nratl were mixed with sample buffer containing lOOmM 
Tris-HCl (pH 6.8), 4% (w/v) SDS, 20% (w/v) glycerol, 0.01% (w/v) 



bromophenol blue, and 200 mM (3-mercaptoethanol. The mixture 
was allowed to incubate at 65 °C for lOmin and SDS-PAGE was run 
using 5-20% gradient polyacrylamide gels (ATTO, Tokyo, Japan). 
After electrophoresis, the gel was transferred to polyvinylidene dif- 
luoride membrane with a semidry blotting system for Western blot 
analysis. The membrane was treated with the primary antibody anti- 
Nratl or anti-H + ATPase (a plasma membrane protein as an internal 
standard) diluted at 1:500 for lh. Then the membrane was washed 
and incubated with anti-rabbit IgG HRP conjugate (1:10 000 dilu- 
tion, Promega) for 1 h at room temperature. The ECL Plus Western 
Blotting Detection System (GE Healthcare) was used for detection 
of signals via chemiluminescence. 



Sequencing of promoter and open reading frame of Nratl 

To compare the promoter sequences, nine rice cultivars were used: 
Kasalath and three cultivars (Miryang-23, PI3 12777, IR36) of the 
Indica subspecies and five cultivars (Koshihikari, Nipponbare, 
Nagoyashiro, Mack Kheua and Mu Bang) of the Japonic subspe- 
cies. The 2.1-kb region upstream of the initiation codon of Nratl 
was amplified by PCR from their genomic DNA using primers 
S'-GGTACCAACACGTCTGACGCTTGTT-S' and 5'-CTCGAG 
ATTCTATGTTGCTAATGCACCTTGT-3' (forward and reverse, 
respectively). To clone the full-length ORF sequence of Nratl, total 
RNA was extracted from rice roots using a RNeasy Plant Mini Kit 
and then converted to cDNA using the protocol supplied by the 
manufacturers of Superscript II. The full-length ORF was ampli- 
fied by PCR using the primers used for amplification of Nratl. The 
fragments of the amplified promoters and cDNAs were cloned into 
the vector pGEM-T (Promega) and the sequences were confirmed 
using a sequence analyser (ABI PRISM 310 Genetic Analyzer, 
Applied Biosciences). 



Genotypic difference in Nratl expression andAI tolerance 
To determine the expression level of Nratl among nine rice culti- 
vars as described above, 5-d-old seedlings were exposed to a 0.5 mM 
CaCl 2 containing 50 uM Al for 8 h. After treatment, the roots from 
each cultivar were sampled for RNA extraction and used for quanti- 
tative analysis of gene expression, as described above. 

To compare Al tolerance, 5-d-old seedlings of nine rice cultivars, 
as described above, were exposed to 0.5 mM CaCl 2 containing 0 or 
50 uM Al for 24 h. The root length was measured with a ruler before 
and after treatment. Ten roots were measured for each treatment. 



Results 

Nratl is responsible for Al QTL on chromosome 2 

Three QTLs controlling Al tolerance were previously detected 
on chromosomes 1, 2, and 6 by using 183 backcross inbred 
lines derived from Koshihikari (Al tolerant) and Kasalath (Al 
sensitive) (Ma et al, 2002). A recently identified Al-tolerance 
gene, Nratl, is located to the QTL region on chromosome 
2. To examine whether Nratl is responsible for the QTL, this 
work obtained two chromosome segment substitution lines, 
which contained the Nratl allele from Kasalath (SL204) 
or not (SL205 as a negative control) in Koshihikari back- 
ground (Supplementary Fig. SI available at JXB online). 
First, Al tolerance to different Al concentrations was com- 
pared among Koshihikari, Kasalath, SL204, and SL205. In 
the absence of Al, the four lines showed similar root elonga- 
tion (Fig. 1A). However, in the presence of 30 and 50 uM Al, 
relative root elongation was significantly lower in Kasalath 
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Fig. 1. Physiological characterization of two chromosome segment substitution lines. (A, B) Root elongation in 5-day-old seedlings of Koshihikari, 
Kasalath, SL204, and SL205 over 24 h; without Al (A; 0.5 mM CaCI 2 , pH 4.5) and with Al (B; 0.5 mM CaCI 2 , pH 4.5, with 0, 30, or 50 nM Al); root length 
was measured before and after treatment and elongation relative to root elongation without Al was calculated; data are mean±SD of six biological 
replicates. (C, D) Al concentration in cell sap (C) and cell walls (D) of root tips; roots of Koshihikari, Kasalath, SL204, and SL205 were exposed to 0.5 mM 
CaCI 2 (pH 4.5) containing 50 |xM Al for 8h; root tips (0-1 cm) were excised and Al concentration was determined by ICP-MS; data are mean±SD of three 
biological replicates. Different letters above the bars indicate significant differences (P<0.05, Tukey's test). 



and SL204 than in Koshihikari and SL205 (Fig. IB). Relative 
root elongation of SL204 was higher than that of Kasalath, 
indicating that other QTLs also contribute to Al tolerance in 
Koshihikari. This is consistent with this study group's previ- 
ous results (Ma et al, 2002). 

Nratl is an Al transporter localized to the plasma mem- 
brane of the root cells and knockout of Nratl results in 
decreased Al uptake (Xia et al, 2010). If Nratl is responsible 
for the QTL on chromosome 2, SL204 should show a phe- 
notype similar to the knockout line. Therefore, Al concen- 
trations in the root cell sap were compared. At the root tip 
(0-1 cm), Al concentration in root cell sap was significantly 
lower in Kasalath and SL204 than in Koshihikari and SL205 
(Fig. 1C) and Al content in cell walls was higher in Kasalath 
and SL204 than in Koshihikari and SL205 (Fig. ID). 

Furthermore, this work compared the expression level of Nratl 
in the root tips of Koshihikari, Kasalath, SL204, and SL205. In 
the absence of Al, the expression level of Nratl was almost dou- 
ble in Koshihikari than in Kasalath at root tips (Supplementary 
Fig. S2 available at JXB online). In the presence of Al, Nratl 
was induced in both Kasalath and Koshihikari, but Kasalath 



showed expression of Nratl lower than Koshihikari (Fig. 2A). 
Substitution of the Nratl allele from Koshihikari with that from 
Kasalath resulted in decreased expression of Nratl to a similar 
level of Kasalath (SL204, Fig. 2A), whereas substitution of other 
regions did not affect the expression of Nratl (SL205, Fig. 2A). 

Expression analysis of ART1 

A previous study has shown that the expression level of 
Nratl is regulated by ART1 (Xia et al, 2010). To investigate 
whether differential expression of Nratl is caused by varia- 
tions in ART1 expression between Kasalath and Koshihikari, 
this work compared the expression of ART1 in two varieties. 
There was no difference in the expression level of ART1 in the 
root tips between Kasalath and Koshihikari (Fig. 2B). 

Comparison of the promoter and coding region 
sequences of Nratl between Koshihikari and Kasalath 

To understand the mechanisms underlying different expres- 
sion of Nratl in two varieties, this work compared the 
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Fig. 2. Expression of Nratl (A) and ART1 (B) in roots of two chromosome segment substitution lines; 7-d-old rice seedlings were exposed to 50 |xM Al 
for 8 h; RNA was extracted from root tips (0-1 cm) and HistoneH3 was used as an internal standard. Data are mean±SD of three biological replicates. 
Different letters above the bars indicate significant differences (P<0.05, Tukey's test). 



promoter sequence of Nratl up to 2.1 kb upstream from the 
translation initiation site. Five nucleotide substitutions and 
a 5-bp insertion were found in the promoter region of Nratl 
of Kasalath (Supplementary Fig. S3 available at JXB online). 
However, the number of ris-acting elements of ART1 was the 
same in the promoters between two varieties (Supplementary 
Fig. S3 available at JXB online). 

This work also compared the sequence of the Nratl cod- 
ing region between Koshihikari and Kasalath. There were 
five nucleotide differences, resulting in four amino acid 
changes at positions E120K, V326I, T500M, and V515A 
(Supplementary Fig. S4 available at JXB online). Predication 
with SOSUI showed that these amino acid changes did not 
affect the topology of Nratl, which contains 10 transmem- 
brane domains (Supplementary Fig. S5 available at JXB 
online). 

Localization of Nratl protein 

Nratl has been reported to be localized to all cells of root 
tips except epidermal cells in a Japonica variety, Nipponbare 
(Xia et al, 2010). To examine whether the differences in the 
promoter and coding region of Nratl affected localization in 
tissues, this work performed immunostaining using an anti- 
body specifically against Nratl. In the root tip (2mm from 
the root tip), Nratl was localized to all root cells except the 
epidermal cells in both Koshihikari and Kasalath (Fig. 3). 
Furthermore, the signal was enhanced in the roots exposed 
to Al in both varieties. Compared with Kasalath, the signal 
in the Al-exposed roots of Koshihikari was stronger (Fig. 3). 

Transport activity of Nratl for Al in yeast 

To test whether amino acid substitutions in the Nratl protein 
affect the transport activity for Al, this work expressed differ- 
ent Nratl alleles in yeast and compared Al uptake. Expression 
of Nratl alleles from both Kasalath and Koshihikari did not 
affect yeast growth in the absence of Al (Fig. 4A). However, 
in the presence of Al, expression of Nratl from either 
Koshihikari or Kasalath increased the sensitivity of yeast to 



Al toxicity (Fig. 4A), indicating that Nratl from Kasalath 
also has transport activity for Al. 

To quantify Al uptake, Al content in yeast expressing dif- 
ferent Nratl alleles was determined. The Al uptake ability 
of Nratl from Koshihikari was 34% higher than that from 
Kasalath (Fig. 4B). Since this difference might result from 
different expression levels of Nratl proteins in yeast, protein 
levels were determined by Western blot with an antibody spe- 
cific to Nratl. The results showed that yeast expressing Nratl 
from Koshihikari and Kasalath produced similar levels of 
Nratl (Supplementary Fig. S6 available at JXB online). 

Since Nratl proteins from Koshihikari and Kasalath dif- 
fer in four amino acids, this work performed a site-directed 
mutation analysis to examine which amino acid is impor- 
tant for Al transport activity. It was found that a mutation at 
the position of 326 from V to I seemed to slightly reduce Al 
transport activity (Fig. 5), whereas other mutations (E120K, 
T500M, and V515A) did not affect Al transport activity. 

Genotypic differences in Nratl expression and 
promoter and ORF sequences 

To further investigate whether Nratl is responsible for genotypic 
difference in Al tolerance in other varieties, the expression levels 
and Al tolerance were compared in nine varieties differing in Al 
tolerance. A good correlation (r=0.82) was found between Nratl 
expression level and Al tolerance (relative root elongation), irre- 
spective of the subspecies (Fig. 6). 

This work also compared the promoter (2.1kb) sequence in nine 
varieties. However, there was no consistent trend between promoter 
sequence and Nratl expression (Supplementary Fig. S3 available at 
JXB online). Also, no trend was found between amino acid sequence 
and Al tolerance in nine varieties (Supplementary Fig. S4 available 
at JXB online). 



Discussion 

Nratl belongs to the Nramp (natural resistance-associated 
macrophage protein) family but, unlike other members, it 
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Fig. 3. Localization of Nratl in roots. Immunostaining with an anti-Nrat1 antibody was performed in 2 mm from the root apex of Kasalath (A, B) and 
Koshihikari (C, D) in roots exposed for 12h to 0.5 mM CaCI 2 without Al (A, C) or 0.5 mM CaCI 2 containing 30 nM Al (B, D). Bars, 100 |im. 
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Fig. 4. Transport activity of Nratl for aluminium in yeast. (A) Effect of Nratl expression on Al tolerance; S. cerevisiae BY4741 carrying empty vector pYES2, 
Nratl (Kasalath), or Nratl (Koshihikari) were spotted on LPM without uracil medium (from left to right: 1 0 [i\ cell suspension with OD 0.2, 0.02, 0.002, and 
0.0002) with or without AICI 3 at different dilutions and incubated at 30 °C for 3 d. (B) Transport activity of Nratl (Kasalath) and Nratl (Koshihikari) for Al 3+ ; 
yeast cells expressing different Nratl alleles were exposed for 4 h to 50 nM AICI 3 (pH 4.2); Al in the yeast was determined by ICP-MS after digestion with 
2M HCI; data are mean±SD of three biological replicates. Different letters above the bars indicate significant differences (P<0.05, Tukey's test). 



transports trivalent Al (Xia et al, 2010). Knockout of Nratl 
results in decreased Al tolerance of rice (cv. Nipponbare, 
a Japonica variety) and decreased Al uptake (root cell sap 
concentration) but increased Al binding to the cell wall (Xia 
et al, 2010). Therefore, Nratl is involved in Al tolerance 
in rice by transporting Al into the root cells prior to final 
detoxification in the vacuoles, which is mediated by a tono- 
plast-localized transporter, OsALSl (Xia et al, 2010; Huang 
et al, 2012). 



The location of Nratl is within the region of a QTL for Al 
tolerance on chromosome 2 (Ma et al, 2002; Yamaji et al, 
2009). Furthermore, another recent study with genome- 
wide association also showed that variations in Nratl are 
probably responsible for the genotypic variation in Al tol- 
erance in the aus subpopulation (Famoso et al, 2011). The 
present study confirmed this possibility by using chromo- 
some segment substitution lines. When the Nratl region of 
Koshihikari (Al-tolerant) was substituted by that of Kasalath 



Genotypic difference in Nratl expression in rice | 4303 



10 

iz 

8 



c 

C/) 
CD 
CD 

c 2 




A # ^ ^ 



,<b N 



Fig. 5. Site-directed mutagenesis analysis. Yeast cells expressing Nratl 
(Kasalath), Nratl (Koshihikari), or mutated Nratl (E120K, V326I, T500M, 
or V515Afrom Koshihikari) were exposed for 4h to 50 \M AICI 3 (pH 4.2). 
Al concentration in cells was determined by ICP-MS after digestion with 
2M HCI. Data are mean±SD of three biological replicates. Different letters 
above the bars indicate significant differences (P<0.05, Tukey's test). 



I 



S 
cd 



.2 2 

w 

£ 

CL 
X 
CD 

CD 1 
> 

CD 

or 



R 2 = 0.66711 



4>T 



l 5- 



20 40 60 80 

Relative root elongation (%) 



100 



Fig. 6. Correlation between Nratl expression and Al tolerance in nine 
varieties. Seedlings of nine rice cultivars were exposed for 8h to 50 \M Al. 
Expression of Nratl in the roots was determined by quantitative real-time 
PCR. HistoneH3 was used as an internal standard; expression relative to 
Kasalath is shown; data are mean±SD of three biological replicates. Al 
tolerance was evaluated by exposing the seedlings to 0 or 50 |xM Al for 
24 h and root length was measured with a ruler before and after treatment; 
data are mean±SD of 1 0 roots. 



(Al-sensitive), the expression of Nratl was reduced (Fig. 2A) 
and Al tolerance was also decreased (Fig. IB). Furthermore, 
the Al concentration in the root cell sap was decreased 
(Fig. 1C) but that in the cell wall was increased (Fig. ID). 
These phenotypes are similar to a Nratl mutant, although 
the extent of these phenotypic changes is different (Xia et al., 
2010). Furthermore, there was a good correlation between 
Nratl expression and Al tolerance in other varieties differing 
in Al tolerance (Fig. 6). Taken together, these results support 
that Nratl is responsible for the QTL detected on chromosome 
2. In addition, SL204 was more tolerant to Al compared with 
Kasalath (Fig. IB), suggesting that other QTLs controlling 
Al tolerance exist in Koshihikari. This is in agreement with 



previous results that have shown that there are two QTLs for 
Al tolerance on chromosomes 1 and 2 detected in Koshihikari 
(Ma et al, 2002). 

Variations in Nratl may result from protein mutation and 
different gene expression. The localization of Nratl protein 
to the root was similar in Koshihikari and Kasalath (Fig. 3), 
but there are four amino acid differences in Nratl between 
Koshihikari and Kasalath, at positions E120K, V326I, 
T500M, and V515A (Supplementary Fig. S4 available at JXB 
online). These changes do not seem to affect the membrane 
topology of Nratl (Supplementary Fig. S5 available at JXB 
online). When the Nratl allele from Kasalath was expressed in 
yeast, it also showed transport activity for Al (Fig. 4B), indi- 
cating that the loss of function does not result from variation 
in Nratl. Quantitative analysis showed that Al uptake medi- 
ated by the Nratl allele from Koshihikari was slightly higher 
than that from Kasalath (Fig. 4B). Furthermore, an amino 
acid change at the position of 326 from V to I seems impor- 
tant to Al transport activity (Fig. 5). Genome-wide associa- 
tion analysis (Famoso et al, 2011) suggested that an amino 
acid change from valine to alanine (but should be threonine 
to methionine, amino acid 500 in the paper) is responsible for 
natural variation in Nratl because this amino acid site was 
predicted to be involved in PKA-type AGC kinase phospho- 
rylation. However, site-directed mutation analysis showed that 
this mutation did not affect the transport activity for Al in yeast 
(Fig. 5). Overall, alterations in protein sequence are unlikely 
to be a major factor responsible for genotypic difference in 
Al tolerance because these alterations were not found in other 
Al-sensitive varieties (Supplementary Fig. S4 available at JXB 
online). However, the possibility could not be excluded that 
all four mutations are required for a difference in Al transport 
activity, which remain to be examined in future. 

By contrast, the expression level of Nratl is responsible for 
the genotypic differences in Al tolerance. This is supported 
by a good correlation between Nratl expression level and 
Al tolerance in different rice varieties (Fig. 6). The expres- 
sion of Nratl is regulated by ART1 (Yamaji et al, 2009; Xia 
et al, 2010), therefore, different expression of Nratl may be 
caused by variation of ART1. However, there was no differ- 
ence in ART1 expression between Koshihikari and Kasalath 
(Fig. 2B). The numbers of c/s-acting elements recognized by 
ART1 in the Nratl promoter region were also similar between 
two varieties (Supplementary Fig. S3 available at JXB online). 
These results indicate that ART1 is unlikely to be a factor 
responsible for the genotypic difference in Nratl expression. 

Alterations in the promoter region have been revealed to be 
the mechanisms for regulating the expression of Al-tolerance 
genes in several studies (Delhaize et al, 2012; Ma et al , 2014). 
For example, the high expression level of ALMT1 is associ- 
ated with tandem repeated elements in the promoter region 
in most Al-tolerant lines of wheat (Sasaki et al, 2006; Ryan 
et al. , 20 1 0). In sorghum, tourist-like miniature inverted-repeat 
transposable elements in the promoter region of SbMATE 
were suggested to be involved in regulating the expression of 
this gene (Magalhaes et al, 2007). In barley, a 1-kb insertion 
(CACTA-like transposon) in the 5'-untranscribed region of 
HvAACTl enhances its expression in Al-tolerant accessions 
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(Fujii et al, 2012). Recently, the higher expression level of 
TaMATElB in several Brazilian wheat lines was found to be 
associated with the presence of a Sukkula-like transposable 
element (11.1 -kb) in its promoter (Tovkach et al. , 20 1 3). More 
recently, in Yorkshire fog (Holcus lanatus), higher HIALMT1 
expression was achieved by an increase in the number of 
cK-acting elements for transcription factor HIART1 in the 
promoter region (Chen et al, 2013). Maron et al. (2013) 
found that Al-tolerant cultivars of maize have three copies 
of ZmMATEl in the genome, which are identical and part 
of a tandem triplication. Therefore, this work compared the 
promoter sequences of Nratl in Al-sensitive and Al-tolerant 
rice varieties. However, no clear differences were found in the 
promoter region up to 2.1kb (Supplementary Fig. S3 avail- 
able at JXB online). There are two possibilities for this result: 
either the regulatory factor is present further upstream of the 
promoter region or the expression of Nratl may be regulated 
differently. Further work is required to understand the mech- 
anism underlying the Nratl expression. 

In conclusion, these results show that differential expression 
of Nratl is mainly responsible for the Al-tolerance QTL located 
on chromosome 2 in rice. Enhancement of Nratl expression in 
Al-sensitive varieties may increase their Al tolerance. 
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